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We investigate an analytic model to compute nonlinear power spectrum of dark matter, galaxies 
and their cross-correlation. The model is based on Press-Schechter halos, which cluster and have 
realistic dark matter profiles. The total power spectrum is a sum of two contributions, one from 
correlations betwen the halos and one from correlations within the same halo. We show that such a 
model can give dark matter power spectra which match well with the results of N-body simulations, 
provided that concentration parameter decreases with the halo mass. 

Galaxy power spectrum differs from dark matter power spectrum because pair weighted number 
of galaxies does not scale with the halo mass and because most halos harbor a central galaxy. If 
the pair weighted number of galaxies increases less rapidly than the halo mass, as predicted by 
theoretical models and observed in clusters, then the resulting power spectrum becomes a power 
law with the slope closed to the observed over several orders of magnitude in scale. Such a model 
also predicts later onset of nonlinear clustering compared to the dark matter, which is needed to 
reconcile the CDM models with the data. Generic prediction of this model is that bias is scale 
dependent and nonmonotonic. This is particularly important for red or elliptical galaxies, which 
are preferentially found in larger mass halos and for which bias in power spectrum may be scale 
dependent even on large scales. 

Our predictions for galaxy-dark matter correlations, which can be observed through the galaxy- 
galaxy lensing, show that these cannot be interpreted simply as an average halo profile of a typical 
galaxy, because different halo masses dominate at different scales and because larger halos host 
more than one galaxy. We compute predictions for the cross-correlation coefficient as a function of 
scale and discuss the prospects of using cross-correlations in combination with galaxy clustering to 
determine the dark matter power spectrum. 



I. INTRODUCTION 



Correlations in dark matter contain a wealth of information about cosmological parameters. Their power spectrum 
is sensitive to parameters such as matter density, Hubble constant, primordial power spectrum slope and amplitude, 
massive neutrinos, baryon density etc. Determining the linear power spectrum of dark matter is one of the main 
goals of modern cosmology. There are several complications that prevent us at present from reaching this goal. First, 
on small scales the linear power spectrum is modified by nonlinear evolution which enhances its amplitude over the 
linear spectrum. It is important to understand this process, so that one can predict the relation between the two. 
This is necessary both to reconstruct the linear spectrum from a measured nonlinear one and to verify whether there 
are other mechanisms besides gravity that modify the clustering of dark matter on small scales. Examples of such 
are baryonic feedback effects on dark matter pj or nongravitational interactions between dark matter particles 0. 
Second, it is difficult to observe correlations in dark matter directly. Direct tracers such as peculiar velocity flows or 
weak lensing still suffer from low statistics and poorly understood systematics. Instead it is much easier to observe 
correlations between galaxies @] or correlations between galaxies and dark matter While these are related to the 
dark matter correlations, the relation may not be simple. The goal of this paper is to address both issues with a model 
that is simple enough to allow analytic calculations without the use of N-body simulations, yet sufficiently accurate 
to be useful for predicting galaxy and dark matter power spectrum. 

Our approach to dark matter clustering is based on the Press & Schechter model j(| . In this picture at any given 
time all the matter in the universe is divided into virialized halos. These halos are correlated and have some internal 
density profile, which can be a function of halo mass. By specifying the halo mass function, their clustering strength 
and their halo profile we can determine the dark matter correlation function. The formalism for correlations inside 
halos has been developed by Q and applied to power law halos ||. We generalize this approach by including the 
correlations between halos and by using more realistic non-power law halo profiles whose shape may depends on the 
halo mass Q. We show in this paper that such a gen eralized model can provide very good agreement with results of 
numerical simulations over a wide range of scales |10(| . 

The central question in extracting dark matter power spectrum from that of the galaxies is how well galaxies 
trace dark matter, the issue of bias. This has been addressed theoretically both with hydrodynamic |ll|] and semi- 
analytic methods jlgJOj ]. The fact that the galaxy correlation function is a power law over several decades in scale, 
while power spectra in CDM models do not show such behaviour, already indicates that the bias is scale dependent. 
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Moreover, galaxies come in different types and observational data show that they can be biased relative to one another 
||l4|| . In our modelling of galaxy correlations we introduce two new functions, the mean number and the mean pair 
weighted number of galaxies inside the halo as a function of the halo mass. The importance of these has recently been 
emphasized in the context of pairwise velocity measurements [ p^|Jl6| and galaxy clustering [12| . These play a key role 
in understanding the relation between galaxy and dark matter clustering. We explore the predictions for different 
choices of these relations and compare them to the results of semi-analytic models. 

Galaxy-dark matter correlations can provide additional information on the clustering of galaxies and dark matter 
and the relation between them. Such correlations have been observed through gravitational lensing effects, for example 
using galaxy-galaxy lensing or correlations between foreground and background populations . Such measurements 
are often interpreted either in terms of an averaged density profile of a halo Q or in terms of a constant bias model 
p8| . We discuss the applicability of these models and how they can be generalized to take into account effects such 
as broad range of halo masses and multiple galaxies inside halos. 



II. DARK MATTER POWER SPECTRUM 



The halo model for power spectrum assumes all the matter is in a form of isolated halos with a well defined mass 
M and halo profile p(r,M). The halo profile is defined to be an average over all halos of a given mass and does 
not necessarily assume all halos have the same profile. The mass is determined by the total mass within the virial 
radius r Vl defined to be the radius where the mean density within it is S v - lr times the mean density of the universe. 
Throughout the paper we will use ACDM model with Vl m = 0.3, Ha = 0.7, normalized to <r$ = 0.9 today. For this 
model <5 v j r ~ 340, although we will also use i5 v ; r ~ 200 (the value for Einstein-de Sitter universe) for consistency with 
the results of some of the N-body simulations. The halo profile is spherically averaged and assumed to depend only 
on the mass of the halo. We will model the halo density profile in the form 

P[r) = (r/r s )-«(l + r/r s )3+«- (1) 

This model assumes that the profile shape is universal in units of scale radius r s , while its characteristic density p s 
at r s or concentration c = r v /r s may depend on the halo mass. The halo profile is assumed to go as r~ 3 in the 
outer parts and as r a in the inner parts, with the transition between the two at r s . The outer slope is fixed by 
the results of N-body simulations which generally agree in this regime. An example of such a profile is a = — 1 
P; ^9|j2C[| . Other models have however been proposed with a = —1.5 |^l|,^2| or even a > — 1 [f23fl . In principle a 
could also be a function of mass scale and may steepen towards smaller mass halos with a ~ — 1 for cluster halos 
and a ~ —1.5 for galactic halos |p4|. Similarly, concentration c may depend on the mass and different authors find 
somewhat different dependence |9j j21| , ^5| . We will explore how variations in the profile and concentration affect the 
power spectrum. Instead of r s we use the concentration parameter c = r v /r s as a free parameter. Note that r v is 
related to M via M — 47r/3r^<5 v irp- Similarly we can eliminate p s and describe the halo only in terms of its virial 
mass M and concentration c, because the integral over the halo density profile (equation [T]) must equal the halo mass. 

For a complete description we need the halo mass function dn/dM, describing the number density of halos as a 
function of mass. It can be written as 

^-dM = ^-J{v)dv, (2) 
dM M w ' w 

where p is the mean matter density of the universe. We introduced function f(v), which can be expressed in units in 
which it has a universal form independent of the power spectrum or redshift if written as a function of peak height 

u=[5 c (z)/a(M)] 2 . (3) 

Here 6 C is the value of a spherical overdensity at which it collapses at z (S c — 1.68 for Einstein-de Sitter model) and 
a(M) is the rms fluctuation in spheres that contain on average mass M at initial time, extrapolated using linear 
theory to z. The form proposed by Press & Schechter (PS) || is vf(v) — (is/2ir) 1 / 2 e~ 1 '/ 2 . This has been shown to 
overpredict the halo abundance by a factor of 2 at intermediate masses below nonlinear mass scale M* p^ , ^7| . A 
modified version of this form that fits better the N-body simulations is given by Sheth & Tormen (ST) [|27fl 

uf(y) = A(l + v'-v)v lX l 2 e- v '' 2 , (4) 

where v' — av with a = 0, 707 and p — 0.3 as the best fitted values, which gives vf{v) oc v 2 for small v. PS 
expression corresponds to a = 1, p = giving vf(y) oc v 5 for small v. The constant A is determined by mass 
conservation, requiring that the integral over the mass function times the mass gives the mean density 



2 



Note that we can still apply this equation even if some dark matter is not bound to any halo. In this case the mass 
function has a nonvanishing contribution in the limit M — > 0. 

The correlation function consists of two terms. On large scales the halos are correlated with each other. We assume 
the halo-halo correlation function follows the linear correlation function. Its amplitude depends on the bias for each 
halo. Halos more massive than the nonlinear mass scale M* are more strongly clustered than the matter, while those 
with masses below M* are less strongly clustered than the matter. A simple halo biasing scheme has been given by 
p8| , ^9| and generalized to the ST mass function by 

Since halos are not pointlikc we need to convolve the halo-halo correlation function with the halo profiles of both 
halos to obtain the dark matter correlation function. The expressions simplify significantly in Fourier space, where 
convolution becomes a multiplication with the Fourier transform of the halo profile 

p(k, M) = [ 4irr 2 drp(r, M) S J^jM . (7) 
J kr 

Note that this is normalized so that p(0, M) = M . It is convenient to renormalize it to unity by introducing a new 
variable y(k,M) = p(k,M)/M, so that y(0, M) = 1 and y(k > 0, M) < 1. The mass of the halo rapidly increases 
as r 3 ~ a up to r = r s , but increases only logarithmically between r s and r v if the outer profile is p(r) oc r~ 3 . The 
dominant contribution to the mass therefore comes from radii around r s . For kr s <C 1 we have y m 1. At kr s ~ 1 
there is a transition and y begins to decrease with k, so that for kr s > 1 we have y(k,M) oc (kr s )~( 3+a \ 

Because the expressions simplify significantly in Fourier space we will in the following only describe the power 
spectrum analysis. The halo-halo term is given by the integral over their mass function with the appropriate bias and 
the halo profile transform, 



PZ(k) = flint*) 



f(u)dub(u)y[k,M(u)} 



(8) 



where Pi; n (fc) is the linear power spectrum and M is related to v via equation^ using the relation between <j 2 {M) — 
4ir J P\i n (k)W B.{k)k 2 dk and M = 4irR 3 p/3, where Wn(k) is the Fourier transform of the top hat window with radius 
R. This gives M oc zy 3 /(™+ 3 ), where n is the slope of the linear power spectrum at scale k ~ We can also define 

the nonlinear mass scale M* where v = 1. Note that on galaxy and smaller scales n < —2 and the relation between 
M and v is very steep, M oc v 1 with 7 > 3. The requirement that on large scales (k — > 0, y ~ 1) the power spectrum 
reduces to the linear power spectrum imposes a nontrivial constraint on the bias distribution, 

f{v)dvb(v) = 1. (9) 

This implies that if halos are biased (b > 1) for M > M* at least some of the halos with M < M* must be antibiased 
(6 < 1) to satisfy this constraint. Most of the bias descriptions in the literature satisfy this constraint to within a 
few percent. The halo-halo term follows the linear power spectrum on large scales and drops below it on scales where 
finite extent of the halos become important (ie where y(k,M) < 1). This term is shown in figure [l] and as expected 
is dominant on large scales. 

In addition to the halo-halo correlation term there are also correlations between dark matter particles within the 
same halos. These are expected to dominate on small scales. We denote this the Poisson term, which is given by 



PL(k) = j— 3 J f{v)du^\y[k,M{v)]\\ (10) 

The main difference between this term and the halo-halo term in equation || is that we have an additional mass 
weighting M/p. This makes the dominant contribution to this term to come from the higher mass halos relative to 
the halo-halo term. On large scales (k — > 0, y ~ 1) the Poisson term is independent of k and behaves as white noise. 
It increases with k more rapidly than the halo-halo term, which scales as the linear power spectrum (figure |l|) . The 
Poisson term declines below the white noise on small scales where the effects of the halo profile become important. 
The total power spectrum is the sum of the two contributions, 
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P dm (k) = PZ(k)+PL(k)- (11) 

To complete the calculation we need to model the dependence of c on M. We will parametrize it as 

c=c (M/M*f. (12) 

Typical values for Co are around 10 at the nonlinear mass scale for a = — 1 profile [j9|]25f| and about a third lower 
for a — —1.5 profile j2lj|. Numerical studies also show that the concentration decreases slowly with the halo mass, 
making (3 negative. 

Figure [l] shows the individual contributions and the sum in comparison to the linear power spectrum and the 
nonlinear prediction from J3l| (PD). In top of the figure we used a = —1.5 and c(M) — 6(M/M*) -0 ' 15 . The latter fits 
the concentration mass dependence given in . Note that for consistency with M] we use 5 V i T = 200 in this case as 
opposed to <5 v i r = 340. In bottom of the figure we used the ST mass function and a = — 1 with c(M) = 10(M/M*)~ - 2 , 
which is somewhat steeper than numerical studies predict [ p5[ as discussed below. The agreement in both cases is 
quite remarkable given the simple nature of the model. It correctly predicts the transition between the linear and 
nonlinear power spectrum, as well as reproduces well the slope at higher values of k. This shows that given a suitable 
choice of c{M) both models can reproduce the nonlinear power spectrum. Conversely, the slope of the power spectrum 
at high k is not directly determined by the inner slope of dark matter profiles, at least if the inner profiles are shallower 
than a = —1.5. 
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FIG. 1. Comparison between the power spectrum predicted with our model and the PD nonlinear power spectrum for 
AC DM model. Also shown are the linear power spectrum and the two individual contributions, P F and p hh . Top is the 
a — —1.5 profile, bottom is a — —1. Other parameters are given in the text. 

In the case of a = — 1 profile the best fitted value for cq — 10 agrees well with ^,||, while (3 ~ —0.2 is somewhat 
lower than f3 — —0.07 J^] and = —0.13 p5| . If one adopts such shallow dependence of c(M) with f3 ~ —0.1 then 
for k > 10/iMpc -1 the predictions of the model are systematically below the PD model. Before concluding that 
this is caused by the galactic halos not being sufficiently compact we must investigate the possibility that the mass 
function is underestimated at small masses. Replacing ST with PS does not significantly affect the results. However, 
both PS and ST assume that each mass element belongs to only one halo, counting only the isolated halos. This 
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is certainly a valid description on large scales, where the total halo mass determines the white noise amplitude of 
the power spectrum. On small scales the dumpiness caused by subhalos within the halos may become important. 
Recent numerical simulations have in fact shown that most of the small halos that merge into larger ones are not 
immediately destroyed, but stay around for some time until they are finally merged on the dynamical friction time 
scale |p0| , pl| ) p3| . In such a case a given mass particle can be part of more than one halo at any given time. Because 
on very small scales the correlation function is dominated by the small halos it may make a difference whether the 
mass is smoothly distributed within the halos or some fraction of it is in the subhalos. However, the contribution to 
the total mass of the halo coming from the subhalos is below 10% ^3|,^(J. Recently, the mass function for subhalos 
from high resolution simulations was determined and it was shown that it is an order of magnitude below the one for 
isolated halos ]32j . One may conclude therefore that subhalos do not affect the mass function significantly and cannot 
resurrect (3 > —0.13, a = — 1 model. 

Steepening the halo profile or changing concentration can both increase the power spectrum to agree better with 
N-body simulations (figure [l]). This is because to increase the power on small scales one has to increase the amount 
of mass contained within a given radius. This can be achieved either by making the inner profile steeper than a = — 1 
or making the concentration parameter larger towards the smaller mass halos. The change in the slope would support 
the results in pl| , p2|] , where the universal profile has the inner slope close to a ~ —1.5, or in pi| , where the profile is 
not universal and steepens as the halo mass is decreased, so that the inner slope changes from a ~ — 1 on the cluster 
scales to a ~ —1.5 on the galactic scales. If the inner slope of the halo profile is a ~ — 1 the concentration has a 
stronger mass dependence than in PJ25|, although the discrepancy is not large. As shown in figure [j] both models can 
fit the nonlinear power spectrum on small scales remarkably well. 
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FIG. 2. Contribution to the P P (k) from different halo mass intervals for the two models in figure |l|. Short dashed 
lines from left to right are M > 1O 14 /i _1 M , 1O 14 /i _1 M > M > IQ^H^Mq, 1O 13 /i _1 M > M > lQ 12 ^ 1 Mq and 
1O 12 /i _1 M > M > 10 11 /i^Mq. Solid line is the total P(k), dotted the correlated term P hh (k). 



Further insight into the relation between the halos and the dark matter power spectrum can be obtained by 
investigating the contribution to the power spectrum from different mass intervals. This is shown in figure || for 
the Poisson term, using the two models from figure l] On large scales the Poisson term is dominated by very 
massive clusters with M > 10 14 /i _1 Mq. These halos dominate the nonlinear clustering on scales around and below 
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k < lhMpc 1 . On smaller scales the contribution from large clusters is suppressed because y(k, M) begins to decrease 
from unity at k ~ r~ x oc c(A/)M -1 / 3 oc M~ ' 5 . This occurs at lower k for the higher mass halos. As a result around 
k ~ lOMvlpc" 1 the halos with 1O 14 /i _1 M0 > M > 1O 13 /i _1 M dominate, while around k ~ 100/iMpc" 1 the halos 
with 10 13 /i _1 Mq > M > 10 12 /i _1 Mq dominate. Note again that the inner slope plays a subdominant role in 
determining the amplitude of the power spectrum. Even if for a steeper slope the power spectrum from a given mass 
interval is decreasing less rapidly (for example for a — —1.5 it is asymptotically flat as opposed to decreasing as fc _1 
for a = — 1), when this becomes important the smaller mass halos have already taken over as a dominant contribution 
to the power spectrum. The nonlinear power spectrum therefore does not reflect the inner slope of the halo profile, 
but rather the halo mass function and the radius at which the mass enclosed within this radius begins to deviate 
significantly from the total halo mass. In both models the halos with M > 10 1:l /i _1 Mq dominate the power spectrum 
for k < 100/iMpc -1 . Any modifications in the linear power spectrum on mass scales below M ~ 10 11 /i _1 Mq |33| ] 
would therefore show up in the dark matter correlation function only on kiloparsec scales and below. 

It is interesting to explore in more detail the quasi-linear regime, where P p (k) ~ const. This approximation is 
valid up to 4irk 3 P(k) ~ 10 or k ~ 1/iMpc -1 . On scales larger than these the power spectrum can be approximated 
as a sum of a line ar p ower spectrum and a constant term, whose amplitude is given as an integral over the mass 
function (equation ^| with y = 1). From figure || one can see that the amplitude of this integral is dominated by the 
massive halos, M > lQ lA h~ x MQ, It is important to emphasize that this amplitude depends only on the integral over 
the power spectrum and not on the details of the power spectrum itself. Even if there are sharp features in the linear 
power spectrum, such as for example baryonic wiggles |34|| , these would not show up as features in the quasi- linear 
power spectrum. Instead, they would be integrated over into a single number, corresponding to the mass weighted 
integral over the mass function (equation [To| ) . This argument is in agreement with the results of N-body simulations 
p5| which indeed show that any baryonic features are erased in the nonlinear regime. This suggests that while the 
PD model breaks down for such spectra, our model could also be applied in such a case. This also applies to the 
spectra with truncated power on small scales ]36|]. We plan to investigate this further in the future. 



III. GALAXY POWER SPECTRUM 



We now apply the above developed model to the galaxies. We assume all the galaxies form in halos, which is a 
reasonable assumption given that only very dense enviroments which have undergone nonlinear collapse allow the 
gas to cool and to form stars. The key new parameters we introduce are the mean number of galaxies per halo as a 
function of halo mass, (N)(M), and the mean pair weighted number of galaxies per halo, {N(N - 1)) 1 /2(M). Just 
as in the case of dark matter these functions are well defined even if the assumption that the statistical properties 
of galaxy population depend only on the halo mass and not on its enviroment is not satisfied pTfl , as long as the 
averaging is performed over all possible enviroments. The resulting power spectrum on small scales where the Poisson 
term dominates is independent of this assumption. On large scales where correlations between the halos are important 
violation of this assumption may lead to a change in the strength of the halo-halo term. 

We furthermore assume that there each halo has a galaxy at its center, while the rest of the galaxies in the halos 
are distributed in the same way as the dark matter, so y(k,M) remains unchanged. This is only the simplest model 
and one could easily generalize it to profiles that differ from the dark matter. Any such complications are important 
on small scales, while on large scales (k < 1/iMpc -1 ) all that is relevant is the total number of galaxies inside the 
halo. The normalization equation ^| becomes 

^W=|, (13) 

where n is the mean density of galaxies in the sample. 
The halo-halo correlation term is given by 



P^{k)=P m (k) 



n S f^)d^b{v)y{k,M) 



(14) 



This should be modified somewhat because the central galaxy does not contribute a y(k,M) term, but this is only 
important on small scales where the halo-halo term is negligible. On large scales where y ~ 1 this term gives the 
constant bias model 

Ofc) = W 2 flin(fc), (15) 

where the mean galaxy bias (b) is given by 
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The Poisson term is given by 



M (N(N-l)) 

Jww J 7 /( } 



y(fc,M)|P. 



(16) 



(17) 



We use the approximation with p = 2 if (N(N — 1)) > 1, because in the limit the number of pairs is large it is 
dominated by the halo galaxies, and p = 1 if (N(N — 1)) < 1, because in the opposite limit the number of pairs in 
this case is dominated by the central galaxy paired with a halo galaxy. Following the usual convention |37| we use 
(N(N — 1)) instead of (N 2 ), since we subtract out the shot noise term arising from the discrete nature of galaxies 
(such a term does not depend on the halo profile y(k, M)). Comparing equations M and |l0| with equations 14 and [l?] 
we see that there is no difference between the two only if (N)/M and {N(N — 1))" 2 /M are independent of M, there 
are many galaxies per halo and the galaxies are distributed as the dark matter within the halo. For such conditions 
the power spectrum of galaxies is identical to the power spectrum of dark matter. 
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FIG. 3. Top figure shows {N(N 
curves) and color Mb < —19.5, Mb 
functions divided by Mh/W 13 M Q . 



M [M D /h] 

l)} 1 / 2 and (N) versus M for galaxies selected by absolute magnitude Mb < —19.5 (upper 
My > 0.8 (lower curves) from semi-analytic models. Bottom figure shows the same 



To test the model above we use semi-analytic models of galaxy formation developed in These models use 

N-body simulations to identify the halos and their progenitors. Gas is assumed to follow dark matter initially so that 
it heats up during the collapse to the virial temperature of the halo. Because of the high density it can efficiently 
cool and subsequently concentrate at the center of the halo. Stars are formed from this cold gas on the dynamical 
time scale. The parametrized star formation efficiency and the stellar population synthesis models are used to assign 
magnitudes in different color bands to the formed galaxies. The small halos with galaxies in them subsequently merge 
into larger halos and exist as individual galaxies until they merge with the central galaxy on the dynamical friction 
time scale. The output of these models is a catalog of halos and their masses. For each halo the output consists 
of a list of galaxies, their positions and luminosities in different bands. From such a catalog one can reconstruct 
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the 3-d distribution of galaxies and dark matter, as well as (N) and {N(N — l)) 1 / 2 averaged over a given range of 
halo masses for any desired galaxy selection criterion. The goal of our comparison is to compare the galaxy power 
spectrum predicted from our model using (N)(M) and (N(N — 1)) 1 / 2 (M) from semi-analytic models to the galaxy 
power spectrum obtained directly from these models. This is a meaningful comparison even if semi- analytic models 
do not correctly describe the nature. If we determine that the model contains all the necessary ingredients to predict 
the galaxy correlations we can then try to obtain these ingredients by other means, either through direct observations 
or better modelling. This can also be applied in the other direction: from observations of galaxy power spectrum (and 
galaxy-dark matter power spectrum discussed in the next section) we can determine the ingredients of our model, 
which must be satisfied by any theoretical galaxy formation model. 

A generic outcome of theoretical models such as these is that the amount of cold gas that can be transformed to 
stars increases as a function of the mass slower than the halo mass itself, because larger halos are hotter and the gas 



takes longer to cool ]ll|,|13|l2|. In such models one would expect (N)/M to decline with M. This is shown in figure 



where (N) and (N(N — l)} 1 ' 2 is plotted versus M for galaxies selected only on the basis on absolute magnitude 
~{Mb < —19.5). Both functions have similar dependence for M > 10 14 /i _1 Mq. When the number of galaxies per halo 
begins to drop below unity the two functions begin to deviate from one another and (N(N — l)) 1 / 2 drops below (N). 
This is because only the halos with two or more galaxies contribute to (N(N — l)) 1 / 2 , while single galaxy halos also 
contribute to (N). However, both functions increase less rapidly than the mass for M > 10 13 h Mq. 




FIG. 4. Comparison between galaxy and dark matter power spectrum predictions for galaxies selected by absolute magnitude 
Mb < —19.5 as in figure ^[ Poisson, halo-halo and combined terms are shown for the two spectra. Also shown is the measured 
power spectrum of galaxies. Note that at low k the Poisson term for the galaxies is lower than that for the dark matter and 
this delays the onset of nonlinear clustering in galaxies. 

Using {N(N — V)) 1 1 2 I'M and (N)/M for Mb < —19.5 from figure || in equations |l3T7 we obtain the galaxy power 
spectrum shown in figure ^. We only show results for a = — 1 model, but the a = —1.5 model gives essentially 
identical results. Also shown is the dark matter power spectrum and its two contributions, as well as the measured 
APM and scaled IRAS galaxy power spectrum compiled in Q . First thing to note is the good agreement between our 
analytical model and the simulations. The agreement is significantly better for this model than for the model where 
there is no central galaxy, which would give a stronger decline in power on small scales. The galaxy power spectrum 
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is almost a perfect power law over several decades in scale, in agreement with observations and in contrast to the 
dark matter power spectrum, whose slope gradually decreases with k. The slope of the galaxy power spectrum is in 
agreement with the observed slope k 3 P(k) oc fc 18 and this slope persists in the analytic model down to kpc scales. 

It is useful to introduce bias b(k), defined as the square root of the ratio between galaxy and dark matter power 
spectrum, 

b(k) = [P gg (fc)/Pdm,dm(fc)] 1/2 . (18) 

The bias b(k) is approximately constant and close to unity on large scales, decreases and becomes less than unity 
between 0.3hMpc _1 < k < 6hMpc _1 and then increases for large k. The bias is therefore scale dependent and 
nonmonotonic, both of which as shown below are generic predictions of this model. On very large scales the power 
spectrum is dominated by the correlations between the halos and the internal structure of halos can be neglected. This 
gives the constant bias on large scales, which for the galaxy type considered here is close to unity. On smaller scales 
the halo Poisson term becomes important both for galaxies and dark matter. However, if (N(N — I)) 1 / 2 /M oc 
with ip < the Poisson term for galaxies is lower than the Poisson term for dark matter in the limit y(k,M) = 1. 
This is because the halo Poisson term is larger if halos are rarer. If ip < the dominant contribution in galaxy power 
spectrum is shifted to lower mass halos, which are more abundant and this reduces the Poisson term relative to dark 
matter. Another important factor that reduces the galaxy Poisson term is that (N) exceeds {N(N — l)) 1 ^ 2 below 
M ~ 10 14 /j,~ 1 M (i). (N)(M) determines the mean density of galaxies n in equation [13]. This suppresses the Poisson 
term in equation |17| even if ip = 0. Suppression of the galaxy Poisson term relative to the dark matter delays the onset 
of nonlinear power in the galaxy power spectrum relative to the dark matter, which is clearly seen in figure |j It gives 
a natural explanation for the position of the inflection point in the observed galaxy power spectrum without the need 
to introduce phenomenological double power law spectra ||. While our model is already in a good agreement with 
the data an even better fit would be achieved with a somewhat smaller Poisson term, which would require ip to be 
even lower or (N) to exceed {N(N — l)} 1 ^ 2 even more. This would further delay the onset of the nonlinear clustering. 




k [hMpc- 1 ] 

FIG. 5. Comparison between galaxy and dark matter power spectrum predicted by our model and the results of N-body 
simulations and semi-analytic models. Predictions for galaxies selected by absolute magnitude Mb < —19.5 and Mb — Mv > 0.8 
are shown. 
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On even smaller scales the halo profile y(k,M) becomes important, since it begins to decrease from unity at a 
scale that corresponds to a typical size of the halo, which is smaller for the lower mass halos. Since the galaxy power 
spectrum is weighted towards lower mass halos relative to dark matter the term y(k,M) begins to be important in 
suppressing the Poisson term at a smaller scale. In addition, if each halo hosts a central galaxy then switching to 
p = 1 for small halos also makes the suppression by the halo profile less important. Smaller suppression of the galaxy 
power spectrum relative to the dark matter results in an increase of bias with k. This argues that the decrease of 
b(k) on intermediate scales and the increase on small scales are generic predictions. The overall result of this is an 
approximate power law of the galaxy power spectrum over several decades. Such a power law arises quite generically 
in a CDM family of models where ip < (or (N) > {N(N — I)) 1 / 2 ) and where each halo hosts a central galaxy. We 
note that the latter is required to preserve the power law behaviour to very small scales. A model where p = 2 for all 
halo masses turns below the power law in the power spectrum at k > 50/iMpc -1 , similar to the dark matter. 

The conclusion above that the bias first declines with k and then rises again applies to a normal galaxy population. 
If one selects red galaxies on the basis of color or ellipticals on the basis of morphology then one may expect a different 
bias dependence, since red or elliptical galaxies are preferentially found in more massive halos, such as groups and 
clusters. Figure [| shows that the galaxies selected by Mb — My > 0.8 and Mb < —19.5 are dominant in halos with 
M > 10 14 /i _1 Mq, while their relative fraction declines rapidly below that. The dependence of (N(N — l)) 1 / 2 with M 
is much steeper in this case so that ip > for M < 10 14 ]i~ 1 Mq. In addition (N) ~ (N(N — l)) 1 ^ 2 across the entire 
range of halo masses, a consequence of the fact that most of the red galaxies are not central galaxies, which in these 
models show recent star formation and are therefore not red. 

Figure || shows the comparison between analytic predictions and results from semi-analytic models jn|. We use 
mass dependence from figure || in equations |l3 - 17 for both Mb < —19.5 and Mb — My > .8 g alaxy selection. Also 
shown are the dark matter power spectrum from the model and from the GIF simulations |38j | which were used for 
semi-analytic models. Qualitatively the agreement is excellent, specially for dark matter and Mb < —19.5 galaxies, 
while for red galaxies semi-analytic models predict a somewhat higher amplitude. Part of the disagreement is caused 
simply by dark matter spectrum not being in agreement with PD (our models are chosen so that they agree with 
PD) pq| . We do not show small scales (k > 20/iMpc -1 ) where limited resolution of N-body simulations prevents 
a meaningful comparison. The remaining discrepancy for red galaxies between 0.5/iMpc -1 < k < 20/iMpc -1 can 
only be explained by them not tracing exactly dark matter distribution in halos with M > 10 13 /i _1 Mq. The red 
galaxies must be more centrally concentrated than dark matter in semi-analytic models in order that their power 
spectrum has a higher amplitude than predicted from our model. This is in agreement with direct analysis of galaxy 
distribution inside halos using the same simulations |p9| , where it was found that red galaxies in AC DM model tend 
to be more centrally concentrated that dark matter. Galaxies that form first end up more towards the center of the 
cluster because the violent relaxation during the merging is incomplete. 

In the case of the red galaxies the bias starts with a value larger than unity on large scales. This is because most of 
the red galaxies are in clusters which are biased relative to the dark matter following equation ||. Bias first rises with 
k and then declines. This is just the opposite from the scale dependence of the normal galaxies and is a consequence 
of i)j > for M < 10 14 /i _1 Mq and (N) ~ {N(N — l)) 1 / 2 . This gives rise to the Poisson term larger for the galaxies 
than for the dark matter on large scales. This conclusion is again independent of the distribution of the galaxies 
inside the halos. This is confirmed in figure |5| where on large scales our model agrees very well with the semi-analytic 
predictions. Because the galaxies are preferentially in larger halos relative to the dark matter y(fc, M) suppression is 
more important and the bias declines on smaller scales. This is seen in the power spectrum from the simulations. In 
our model it begins to rise again on even smaller scales because p switches to unity for M < 10 /i _1 Mq, resulting 
in a smaller suppression by the halo profile. This effect is not seen in the simulations, presumably because of their 
limited resolution. 

It is important to note that bias may never be really constant even on scales above 100/i~ 1 Mpc. For the red sample 
it changes by 30% between k = 0.01/iMpc -1 and k — O.l/iMpc -1 . This is because the Poisson term does not become 
much smaller than the halo- halo term even on very large scales, a consequence of the fact that the slope of Pn n (fc) 
and thus the halo-halo term itself becomes flat and even positive on very large scales (approaching n ~ 1 on very 
large scales). Since even at the turnover of the power spectrum (where n ~ 0) the Poisson term for the red galaxies 
is of the order of 20% of the halo-halo term the bias does not become constant and begins to increase again on scales 
larger than the scale of the turnover. In fact on very large scales (k < 10 _3 /iMpc _1 ) the red galaxy power spectrum 
becomes white noise, although these scales are already approaching the size of the visible universe. It should be noted 
that this description is valid on large scales only for galaxies which do not obey mass and momentum conservation. 
For the dark matter mass and momentum conservation require that the Poisson term vanishes on large scales and 
any spectrum generated by a local process should decrease faster than P(k) oc fc 4 as k — + fl40|| . Galaxies do not obey 
mass and momentum conservation and can have the Poisson contribution, so the qualitative scale dependence of bias 
remains as predicted above. 
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We have concentrated on the power spectrum above because it is the quantity that can be most directly compared 
to the theoretical predictions. The same analysis could however be applied to the correlation function as well. The 
power law dependence of the power spectrum would also result in a power law correlation function, so the conclusions 
would remain unchanged. The main difference in the real space is that the Poisson term is localized to scales smaller 
than the typical halo scale and vanishes on scales above that. In this case bias would be scale dependent up to this 
typical scale (of order few Mpc), but would become scale independent on scales above that. There is no need to model 
the Poisson term on large scales at all. In this sense the real space correlation function offers some advantages over 
the power spectrum, where one must attempt to remove the Poisson term in the power spectrum by modelling it as 
a constant term on large scales. 

Our predictions agree with the results of semi-analytic models, indicating that the here proposed model is sufficient 
to extract the key ingredients to model the galaxy clustering. This means one does not need to rely on N-body 
simulations as long as the ingredients of the model are specified. If one can extract (TV), (N(N — l)) 1 / 2 (7\/) and 
y(k,M) directly from the data one can sidestep the theoretical modelling of this relation and predict the galaxy 
power spectrum directly |4l|. It is in principle possible to obtain such information at least for the massive halos by 
combining dynamical information on galaxy groups and clusters, such as X-ray temperature, velocity dispersion or 
weak lensing mass, with the number of galaxies in these clusters. Existing data such as CNOC survey [ }42[ indeed find 
that (N)/M for galaxies with Mk < —18.5 is systematically lower in massive clusters with a > 1000 km/s than in 
poorer clusters. The current data are sparse, but new large surveys such as SDSS and 2dF will enable one to extract 
such information with a much better statistics. This could allow one to determine within our model the dark matter 
power spectrum from the galaxy power spectrum directly. 

Another direction to obtain (N)/M is to require consistency with other measurements that combine dynamical 
and galaxy information. Galaxy-dark matter correlations discussed in the next section are one possibility. Another 
are pairwise velocity dispersion measurements. If (N)/M declines with M then the pairwise velocity dispersion for 
the galaxies will be lower than for the dark matter p^pl. This is because there will be more pairs of galaxies 
in smaller halos relative to the dark matter. Smaller halos have smaller velocity dispersions and smaller relative 
velocities between the particles. This can exp lain the lower amplitude of pairwise velocity dispersion in the LCRS 
data compared to the N-body simulations |lq ]. The required value of ip ~ —0.1 has indeed the same sign as required 
to reproduce the delayed onset of nonlinear clustering and the power law in galaxy power spectrum. It would be 
interesting to see whether a sing le set of functions (N)(M), (N(N - l)) l l 2 {M) can provide a unified description to 
both galaxy clustering and pairwise velocities within the CDM models. 



IV. DARK MATTER-GALAXY CROSS-CORRELATION 



Dark matter-galaxy cross-correlations are measured whenever a galaxy is cross-correlated with a tracer of the 
dark matter. Examples of this are galaxy-galaxy lensing [ fL7| , where one is measuring correlation between galaxies 
and cosmic shear, and correlations between foreground and background galaxies or quasars Q, where correlations 
(or anti-correlations) are induced by magnification bias of background objects. In both cases one is measuring the 
correlations between the galaxies and dark matter along the line of sight, which can be expressed as a convolution 
over the galaxy-dark matter cross-correlation power spectrum. 

Galaxy-dark matter cross-correlations have been modelled in the past using either a bias model relating them to the 
dark matter or galaxy power spectrum |l^] or using galaxies sitting at the centers of the galactic size halos |j| . In the 
first description assuming galaxy-dark matter cross-correlations measure bias 6(fc)Pdm(fc), which in combination with 
the galaxy power spectrum b 2 (k)Pd m (k) can give both b(k) and P& m (k). Such a model is a reasonable description on 
large scales, but must break down on small scales where galaxies do not trace dark matter and there is no guarantee 
that the scale dependent bias that relates P ga i,dm(fc) and P ga i(fc) can be used to extract Pd m (k). 

Second model describes cross-correlations in terms of galaxies sitting at the centers of their halos and interprets 
the results in terms of the averaged halo profile j^]. There are two potential problems with this approach. First, 
there may be more than one galaxy inside the halo, which is specially important for large halos (figure [?]). Since 
not all galaxies can lie at the halo center this can affect the interpretation of the cross-correlations in terms of the 
halo profile. Second, just as in the case of the dark matter the contribution to the power spectrum comes from a 
range of halo masses and one cannot model the galaxy-dark matter cross-correlation simply as a typical L* galaxy 
halo profile. The strength of the correlations is determined both by the dark matter profile of the halos as well as by 
the halo mass function, so the slope of the correlation function that one is ultimately measuring with galaxy-galaxy 
lensing and foreground-background galaxy correlations need not be directly related to the dark matter profile fl43| . 
Model developed in previous sections may be applied to the dark matter-galaxy cross-correlation power spectrum to 
quantify these issues in more detail. 
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Galaxy-dark matter cross-correlation power spectrum has halo-halo and halo Poisson terms. First term describes 
the correlations between galaxies and dark matter in neighbouring halos and is dominant on large scales. Second term 
includes the correlations between the galaxies and dark matter in the same halo and dominates on small scales. The 
halo-halo term is given by 



P - f f(v)dv^±b(v)y[k,M(v) 



f{v)dvb{u)y%M{v)] 



(19) 



On large scales where this term dominates it reduces to constant bias model, Pg h dra (k) = (b)F\ ln (k). The Poisson term 
in the model where galaxies trace dark matter inside the halos except for the central galaxy sitting at its center is 
given by 



1 



f(v)dv(N)\y(k,M)\P. 



(20) 



Here p = 2 for (N) > 1 and p = 1 for (N) < 1. 
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k [hMpcr 1 ] 

FIG. 6. Galaxy-dark matter cross-correlation power spectrum for the two galaxy types as in figure ^ (dashed), together with 
the dark matter (solid) and galaxy (dotted) power spectrum. Also shown is the cross-correlation coefficient r(k) from the model 
(dash-dotted) and from simulations (long-dashed). 



Figure g shows the results for the cross-correlation power spectrum for the same galaxy selection as in figures y and 
|]. For regular galaxies selected by an absolute magnitude (top panel) the cross-correlation spectrum is similar to the 
galaxy power spectrum. If we define the cross-correlation coefficient as 

r (l,\ - PfS-dmjk) , 

m [P dro (fc)P g (fc)]V 2 < (21) 

then we see from figure ^ that it is approximately unity up to k ~ 1/iMpc -1 and increases for higher k. Note that the 
cross-correlation coefficient is not restricted to \r(k)\ < 1 because we have subtracted out the shot noise term from the 
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galaxy power spectrum following the usual approach |37|. Because on small scales the galaxy and cross-correlation 
spectra are comparable and exceed the dark matter spectrum the cross-correlation coefficient grows to large values 
in this model. Comparison with the semi-analytic results [[l3| again shows very good agreement up to the resolution 
limit of the simulations. Bottom of figure [fj] shows the results for the red galaxies. In this case the cross-correlation 
spectrum falls in between the dark matter and the galaxy spectrum, so that r ~ 1 down to very small scales. This is 
again in agreement with semi-analytic results which show r ~ 1 throughout the entire range of k. 

The main reason for r ^ 1 on small scales is that {N(N - l)) 1 / 2 ^ (N) (figure g). The differ ence between the two 
functions is more significant for the normal than for the red galaxies, which is why the cross-correlation coefficient 
begins to deviate from unity at larger scales for Mb < —19.5 than for Mb — My > 0.8. Because in this regime 
{N(N — l)) 1 / 2 < (N) this leads to r(k) > 1, as seen in figure g It is interesting to note from figure gthat for the red 
galaxies the two functions agree very well even below unity and this leads to r(k) ~ 1 down to very small scales. When 
this happens one can reconstruct the dark matter power spectrum from the galaxy and cross-correlation spectrum 
even if most of the dark matter halos are not directly observed. Unfortunately one cannot extract these two functions 
without first identifying the dark matter halos, so this prediction cannot be directly verified from observational data 
using the galaxy information only. 

Second source of stochasticity is the presence of central galaxy. For those halos where {N(N — l)) 1 / 2 < 1 or (N) < 1 
only one power of y(k, M) is used as opposed to two in the case of the dark matter. This induces some stochasticity 
even if (N) = (N(N — l)) 1 / 2 , because it enhances the galaxy-galaxy and galaxy-dark matter spectrum above the dark 
matter-dark matter spectrum. Another source of stochasticity would be ip ^ 0, which would make correlations at a 
given scale being dominated by different mass range in the case of the dark matter and the galaxies. Calculations 
where only this effect is present give r ~ 1 over a wide range of scales, showing that this cannot be a significant source 
of stochasticity, at least for reasonable values of tp. 
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k [hMpc- 1 ] 

FIG. 7. Contribution to the dark matter galaxy cross- correlation power spectrum from the different halo mass intervals. The 
curves correspond to the same mass intervals as in figure ^. Top is a — —1.5 model, bottom a = — 1. 

Our model predicts that even if the constant bias model is not valid, its generalization r = 1 model may be a 
reasonable approximation at least down to 1 Mpc scales. An example are the red galaxies (bottom of figure |(]), 
which have very strong scale dependent bias, yet r w 1 over a wide range of scales. In this sense determining the 
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dark matter power spectrum from the measurements of galaxy-galaxy spectrum and galaxy-dark matter spectrum 
under the assumption of r — 1 may have a larger range of validity than the constant bias model. This relies on the 
assumption (N) — {N(N — l)} 1 ^ 2 predicted from these models. This prediction can be verified at least for the more 
massive halos directly from observations, for example by using galaxy counts in cluster catalogs to extract (TV) and 
(/V(/V - l)) 1 / 2 . Such an approach would provide an alternative way to determine r(k) directly from the data. 

The model developed here can also be used to clarify the interpretation of the galaxy-dark matter cross-correlation 
in terms of an averaged density profile of a typical galaxy. Figure 7^ shows the contribution to the cross-correlation 
spectrum from the different halo mass intervals, similar to figure |2 for the dark matter. For k < 20/iMpc~ 1 , corre- 
sponding approximately to scales larger than 100h _1 kpc in real space, one cannot interpret the correlations in terms 
of the shape of a single halo profile, but instead as the convolution of these over the halo mass function, multiplied 
with the number of galaxies per halo. Observed correlations on large scales do not necessarily mean that the halo of 
an L* galaxy extends to large distances. Instead, it is more likely that one is observing correlations arising from the 
group and cluster size halos, which exceed the correlations contributed from the galactic size halos on larger scales. 
This cannot be corrected in any simple manner by taking into account the correlation function of the galaxies [ |L7| , 
which attempts to model the presence of other nearby halos. Even if the galaxy correlations vanished one would still 
need to take into account the halo mass function and the fact that different halos dominate on different scales. More 
detailed discussion of these points will be presented elsewhere E| . 

On smaller scales the transition to (N)(M) < 1 implies that y(k,M) suppression is less important because p = 1. 
This is further enhanced by the flattening of (N) below M ~ 10 13 /i _1 Mq as seen in figure |3| In addition, galaxies 
selected on the basis of their absolute magnitude cannot exist in very small halos, so the mass function has a strong 
cutoff below lO 12 /?.^ 1 Mq. Thus on scales with k > 20ft.Mpc~ 1 the galaxy-dark matter cross-correlation may be 
better interpreted in terms of the average profile of 10 12 /i~ 1 A/q < M < 10 13 /i _1 Mq halos. However, this may not 
be a robust prediction since the semi-analytic predictions in figure [| are highly uncertain over this mass range. A 
small change in (N}(M) may lead to a larger influence of the mass function on the power spectrum, making the 
correspondence between the halo profile and the power spectrum less certain. In general one should be cautious in 
interpreting the shape of the galaxy-dark matter correlation function in terms of an averaged dark matter profile. 



V. CONCLUSIONS 



We developed an analytic model for computing the power spectrum of the dark matter, galaxies and their cross- 
correlation based on the Press-Schechter model. In this model all the matter in the universe is divided into virialized 
halos. These halos cluster and have some internal profile. The total power spectrum is the sum of the halo clustering 
term and the halo Poisson term, which accounts for the correlations within the halos. We assume that the halo 
profiles are self-similar regardless of the initial conditions, but with the mass dependent concentration parameter, as 
suggested by high resolution simulations |) 23 2l[] . The model agrees well with the results of N-body simulations for 
the KCDM model. We are able to find a good agreement for inner slopes a = — 1 and a = —1.5, indicating that the 
shape of the nonlinear power spectrum cannot by itself distinguish between the two. 

The model can in principle be applied to any cosmological model, including those with a cutoff in the linear power 
spectrum on small scales or with some features in the power spectrum. While this will be explored in more detail in 
a future paper we wish to emphasize here that the mass function, which is sensitive to the linear power spectrum, has 
a direct effect on the nonlinear power spectrum through the halo abundance, so that not all of the information on the 
linear power is lost in the nonlinear regime. For example, if the linear power spectrum is cut-off on small scales and if 
inner profile a > —1.5 as suggested by the simulations then the correlation function or k 3 P(k) must have a turnover 
on small scales. This differs from the CDM models which predict the nonlinear correlation function to continue to 
grow on small scales. If we wish to eliminate the halos with M < 10 11 /i _1 Mq |33j then this would suppress the power 
on scales below lOkpc (figure ||). This effect therefore becomes significant on scales smaller than those resolved in a 
recent study of such truncated power spectrum models |3(| . 

Our main conclusion regarding the galaxy power spectrum is that a simple model for the dependence of the linear 
and pair weighted number of galaxies inside halo as a function of the halo mass can explain most of the properties of 
the galaxy clustering seen in more complicated models based on the N-body simulations. A power law in the galaxy 
correlation function with slope 1.8 is a generic prediction of the model where the number of galaxies inside the halo 
increases less rapidly with mass than the halo mass itself, mean number of galaxies exceeds pair weighted average 
and there is a central galaxy in each halo. The decline of number of galaxies per unit mass as a function of mass 
is predicted by the galaxy formation models p^ , [l3|,|lT[ | and has been observed in clusters [Q. It is also required 
to explain the pairwise velocity dispersion results"||l5fT For such galaxies bias first decreases below unity, because 
the Poisson term is smaller for them than for dark matter. This naturally explains the later onset of nonlinearity 
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in galaxy power spectrum compared to the dark matter, which reconciles the discrepancy between the data and the 
CDM models Conversely, there is no need to invoke poorly motivated models such as double power law model 
j|. On large scales bias converges to a constant for these galaxies. 

Red or elliptical galaxies, which are more abundant in massive halos, show a different relation: their number inside 
the halos increases on average more rapidly than the halo mass. In this case bias increases with k above the turnover 
in the power spectrum (k ~ O.Ol/iMpc^ 1 ), because their Poisson term is larger than that of dark matter. In fact, the 
Poisson term may be so strong that it may not be negligible compared to the halo clustering term even on very large 
scales and one may not converge to the constant bias model. 

Galaxy-dark matter correlations can also be predicted by this model. In this case one must specify the average 
number of galaxies per halo as a function of halo mass. Here again our model reproduces the main features present 
in the N-body simulations with semi-analytic galaxy formation p3j ]. Galaxy-dark matter cross-correlations can be 
measured with galaxy-galaxy lensing or correlations between foreground and background galaxies and may provide a 
way to break some of the uncertainties present with the galaxy clustering. For example, we have shown that even if the 
constant bias may not be a good approximation, cross-correlation coefficient may nevertheless be close to unity down 
to Mpc scales, which would allow one to extract the dark matter power spectrum from the knowledge of the galaxy 
and cross-corelation spectrum on scales larger than this. The main source of stochasticity (r =/= 1) arises from the 
pair weighted number of galaxies inside the halo differing from the mean number of galaxies and from the (possible) 
existence of central galaxies in the halos. 

We have emphasized that caution must be applied when interpreting the cross-correlations such as galaxy-galaxy 
lensing in terms of an averaged density profile of a halo. As we have shown different halo masses dominate on different 
scales and the correlation function reflects this combined effect of all the halos. For example, correlations at a few 
hundred kpc observed by galaxy-galaxy lensing [ p"7[ are more likely to be caused by group and cluster sized halos at 
r s distances than by galaxy sized halos at r v distances. More detailed work is needed to extract the structure and 
extent of the dark matter halos from such observations. 

Perhaps the most promising direction to explore in the future is to extract the functional dependences that 
parametrize our model directly from the observations. If one can determine the linear and pair weighted number 
of galaxies as a function of halo mass and their distribution inside the halos then one can determine the galaxy power 
spectrum directly within this model. Similarly if one can compare the mean number of galaxies with the pair weighted 
number as a function of halo mass then one can predict the galaxy-dark matter cross-correlation coefficient. This is 
certainly feasible for clusters, which dominate the Poisson term on large scales. Current data are sparse (42), but new 
surveys such as SDSS or 2dF should provide sufficient statistics to make this feasible. This approach would provide 
an independent estimate of the scale dependence of bias and correlation coefficient on large scales. It will also provide 
important constraints that would need to be satisfied by any viable galaxy formation model. 
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